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Abstract 

We show that in SU(2) Yang-Mills theories a simple relation exists between lattice gluon prop- 
agators in Coulomb and Landau gauge and discuss the physical implications of such result. In 
particular, the realization of the Gribov-Zwanziger confinement mechanism in Coulomb gauge, 
linked to dual-superconductivity, would imply that the standard BRST charge must be ill defined 
non perturbatively. As a consequence, the Kugo-Ojima confinement criterion, which relies on 
BRST charge conservation beyond perturbation theory, would not be fulfilled. 
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I. INTRODUCTION 



A. Overview 

Understanding confinement in non-Abelian gauge theories has proved a major challenge 
for theoretical physics. How to enforce Gaufi's law, which ensures gauge invariance and 
vanishing color charge on physical states, Q c \^) P h = 0; how to explain the discrepancy 
between color neutral S-matrix hadronic states and color charged fields appearing in the 
action; and how the IR properties of gauge fixed Green's functions should reflect confinement 
of gluons and quarks and non-physicality of ghost fields have been, from the very beginning, 
closely intertwined questions propaedeutic to the appearance of a linearly rising potential 
and which have spawned extensive research dating back to the '70s. 

In a seminal paper, Kugo and Ojima [l| gave a criterion on the IR behavior of the ghost 
Green's function that ensures S-matrix states being color neutral in covariant gauges with 
the Faddeev-Popov (FP) action and the corresponding Becchi-Rouet-Stora-Tyutin (BRST) 
symmetry. According to their result, the ghost dressing function should diverge as the 4- 
momentum goes to zero. The key assumption underlying their derivation is that a global 
BRST charge can be defined beyond perturbation theory, such that QbrstI^OpIi = 0. Color 
neutrality of all S-matrix states and the disappearance of ghosts from the spectrum is then 
ensured by the so-called quartet mechanism [l,]. Equivalently, if Q c \^} P h = and the 
Kugo-Ojima criterion is not fulfilled, BRST symmetry must be non-perturbatively either 
spontaneously broken or an ill-defined concept. 

A complementary approach was pioneered by Gribov [2] and later expanded by Zwanziger 
{3)]. The basic idea is that the FP mechanism is not sufficient to define the non-Abelian 
partition function beyond perturbation theory. To extract a unique field along the gauge 
orbit the functional integral needs to be restricted, via terms which are in general non-local 
and are often referred to as the "Horizon condition", to the Gribov Region f2, where 
the FP operator is positive definite, or even further to the Fundamental Modular Region 
A, where the gauge fixing functional has only absolute minima. Obviously, any restriction 
on the integration domain implemented other than through the standard FP mechanism 
will in general break BRST symmetry. Although it can be applied also to covariant gauges, 
the physical implications of the Gribov-Zwanziger approach are transparent in Coulomb 
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gauge. The latter has the advantage of projecting, in the Hamiltonian formulation, onto the 
subspace of states satisfying GauB's law [J, |5|, so that no additional assumptions to ensure 
Qc\^) P h = 0, as in the Kugo-Ojima approach, are needed and an explicit construction of 
the gauge invariant Hilbert space as in Ref. &] can be circumvented. The Horizon condition 
implies in Coulomb gauge an IR diverging ghost dressing function and a vanishing static 
transverse gluon propagator ^j. Moreover since the inverse of the ghost dressing function 
can be identified with the dielectric function of the Yang-Mills vacuum [7|, the latter behaves 
like a perfect color dia-electric medium, i.e. a dual-superconductor; a confining potential 
can arise via dual-Meissner effect. 



B. Motivations 



Approximation schemes for the set of Dyson-Schwinger equations (DSE) in QCD had 
been proposed already by Mandelstam in the late '70s |8|, leading to the now abandoned 
hypothesis of infrared slavery. A first self consistent solution for the DSE was obtained by 
Cornwall [9;] in the pinch technique framework [loj|, leading to a massive IR behavior of 
the gluon propagator. In such solution the Kugo-Ojima criterion is not fulfilled and BRST 
symmetry is lost non perturbatively. Only from the late '90s on alternative approaches have 
been developed to find solutions to the tower of DSE satisfying BRST invariance Since 
then, advances in functional methods for the non perturbative analysis of Yang- Mills Green's 



functions in Landau 
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13} and Coulomb jlj], ll5j gauge on one side and the availability 

23] 
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of large lattice simulations with improved gauge fixing techniques 
have revived the interest in the IR analysis of gauge fixed Green's functions. 

As mentioned above, in Landau gauge continuum methods allow in general for two types 
of solutions [13]. When coupled to the requirement that a global BRST charge be defined 
and conserved, as assumed by the Kugo-Ojima confinement criterion for the Landau ghost 
dressing function, they lead in any dimension to a conformal behavior of the gluon propagator 
in the IR, the so called "scaling" solution 



11 



13] . Imposing BRST symmetry to hold beyond 



the perturbative regime assumes implicitly that no extra care needs to be taken to restrict 
the functional integral beyond the FP operator. Rela xing this requirement, however, the 
gluon propagator can acquire a massive IR behavior 9], [ljj 24, 25] while the ghost dressing 



function stays finite, also known as "decoupling" solution in the literature. In principle 



BRST symmetry gets explicitly broken when rendering local the terms in the action enforcing 
the Horizon condition, since this introduces additional ghost fields and induces modified 
residual transformations. Such terms should lead to a massive IR gluon propagator {24]. 
Alternatively, it has been argued in the Jiterature that a vanishing gluon propagators can 
still be obtained in this framework 
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28fl and that dimensional two condensates 
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31J 



are eventually responsible for the massive IR solutions as in Ref. 24 1. 

As for the lattice, all results for the gluon and ghost in Landau gauge obtained in "stan- 
dard" setups clearly favor a IR massive gluon in 2+1 and 3+1 dimensional Yang-Mills theory, 
while agreement with scaling is only found in 1+1 dimensions 20|, |32|. Some works apply- 



ing "non-standard" lattice methods have however pointed to possible discrepancies in such 
picture [33|-|37]. 

Whether one or the other scenario is realized in QCD is the subject of an ongoing con- 
troversy. It should be noticed that in both cases the gluon pro pag ator violates reflection 
positivity and disappears from the physical spectrum anyhow 38]. Moreover, a linearly 
rising potential could arise in both frameworks 39[. Indeed, the issue boils down to whether 
color charge descends from BRST charge or, equivalently, whether in Landau gauge the 
Horizon condition can be neglected in favor of the standard FP mechanism. 



In Cou 
results 



omb gauge we are in a more comfortable position: in 3+1 dimensions both lattice 



2l| and functional methods 
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show a IR vanishing static gluon transverse prop- 
agator well described by Gribov's formula ^J, while the Coulomb ghost dressing function 
clearly diverges at zero momentum 41, (42)]. The Gribov-Zwanziger mechanism is realized 
and a coherent picture emerges, naturally enforcing Q c \^) P h = and with a direct phys- 
ical interpretation in terms of dual Meissner effect 7j. As we will show below, the 2+1 
dimensional Yang-Mills theory behaves analogously. 

As for the 1+1 dimensional case, the theory is topological 43] and no physical gluon 



exists. As a "back of the envelope" calculation shows, the Coulomb gauge propagator 
vanishes identically in momentum space and the corresponding FP operator can be inverted 



exactly, while the Landau propagator only describes pure gauge degrees of 



An IR conformal behavior with conserved BRST charge is thus natural 43 



reedom 44 ] . 



44j. Indeed, 



topological theories are the starting points for the construction of non-perturbative BRST 
invariance in any dimensions. Interestingly enough, lattice results, which implement some 
sort of Horizon condition restricting by construction to Q, and continuum calculations using 



only the FP mechanism seem to coincide in this case. 

In higher dimensions however it is not clear whether the Horizon condition in Landau 



gauge can be reconciled with a vanishing g 
arguments have been brought up in favor 



uon propagator. As mentioned above, different 



26M28] and against 
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it. Conventional 



lattice results, including direct simulations of the Kugo-Ojima function [451 ]. are quite unam- 
biguous in favoring the latter approach, although gauge fixing ambiguities and/or problems 
in the definition of the lattice gauge fields have been invoked to claim possible agreement 



with the former scenario 



33, 
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46], though the massive gluon IR behavior turns out to 



be very robust even when choosing ad-hoc prescriptions for the gauge fixing procedure to 
enhance the ghost propagator [36] . It is also not clear at the moment how a redefinition 



of the gau ge fi elds 



46] would circumvent the well known no-go theorem for a lattice BRST 



symmetry 47], |48j . Even if such approach would work, it would anyhow imply a deep change 
of paradigm for lattice gauge fixing, since on one hand a BRST invariant construction would 
sum over all copies in the different Gribov regions while on the other hand the standard 
gauge fixing procedure on the lattice follows Gribov's spirit, effectively restricting the con- 
figurations to lie within Q. There is of course a residual ambiguity in the choice of the best 
among all local minima found, since the explicit restriction to A would require an infinite 
numerical precision, but on one side going to the continuum limit makes the gauge ambi- 
guity much milder 22], while on the other side thermodynamic arguments suggest that the 
relevant contribution to the functional integral will anyway lie on <9A R dQ, so that gauge 
ambiguities, apart from topological ones, should not play any role. 

Since i) reliable physical input through phenomenological data is still lacking; ii) the 
matter turns out difficult to solve by brute force numerics, i.e. by explicitly going to the 
thermodynamic and/or continuum limit and iii) still waiting for a valid and viable BRST 
construction on the lattice, we follow here an independent approach. In dimensions higher 
than 1+1 a physical, although confined, gluon obviously exists and one should be able to 
read its degrees of freedom out of both the transverse Landau and Coulomb propagators: 
a relation between the two and eventually to the physical gluon degrees of freedom should 
exist. Moreover if Gribov's prescription to restrict the functional integral is the physically 
sound one in Coulomb gauge, since it seems to agree with dual superconductivity, one would 
expect the Gribov mass appearing there {2], Q to be linked to some physical quantity, i.e 
it should effectively be gauge invariant. This mass should show up in Landau, interpolating 
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49[ and covariant gauges [50| and be related to the effective IR gluon mass [9]. 

In d+1 dimensions the most natural Ansatz to relate the static transverse gluon prop- 
agator in Coulomb gauge and the transverse gluon propagator in Landau gauge is via a 
one to one correspondence between the variables they must be function of, namely fr 2 , the 
natural 0(d) invariant in the static Hamiltonian Coulomb gauge picture, and p 2 , its 0(d+l) 
covariant counterpart in euclidian lattice Landau gauge. We therefore assume that in 2+1 
and 3+1 dimensions by just relating monotonically the <i-and d + 1-momentum scales the 
two propagators will coincide. Based on this simple Ansatz we will show that i) there is a 
one to one correspondence between the lattice results for the IR vanishing static Coulomb 
gluon propagator, realizing Gribov's original Horizon condition j^j and the IR massive Lan- 
dau propagator and ii) that within this correspondence the Coulomb gauge Gribov mass 
exactly coincides with the Landau gauge effective IR gluon mass. Should the uniqueness of 
the massive behavior for lattice gluon propagator in Landau gauge be confirmed, one could 
interprete the loss of BRST symmetry beyond perturbation theory that would descend from 
it as a consequence of the dual- superconductor mechanism in QCD. 

II. RESULTS 

We calculate the transverse gluon propagator in Coulomb and Landau gauge in 2+1 and 
3+1 dimensional SU(2) Yang-Mills theory on configurations generated via MC simulations 
on lattices up to 64 3 and 32 and /3 between 3-12 and 2.15-2.6 respectively, using a heath 
bath plus over-relaxation algorithm. The scale is set re-expressing the lattice spacing in 



units of the string tension 



3, 3, 



which we fix to a = (440 MeV) 2 also in 2+1 dimensions. 



A. Coulomb propagator 

In Coulomb gauge renormalization issues on the lattice have only been recently clarified 



in Ref. 
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52]. We extend here our procedure to 2+1 dimensions and improve our gauge 
fixing algorithm. 

For each fixed /3 we define the gluon propagator D^(jp) as the Fourier transform of the 
gluon two-point function: 

2 . /7rJfe p 



D$(p) = (At(k)A%-k)) rr = 5 ab 5 l3 D p (p) , p„ =p(*,) = - sin ( ) . (1) 



Here A a ^{k) is the Fourier transform of the lattice gauge potential A a (x + jx/2), defined as: 

A,{x + fi/2) = ^{U,{x)~U,{x)^) , (2) 

a is the lattice spacing, p = (p,Po) denotes the four-momentum and G (— L/2, +L/2] 
are the inte ger- valued lattice momenta. We always select spatial indices k satisfying a 



cylinder cut 
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41] to minimize violations of rotational invariance. Time-like momenta po 
are unconstrained. 

I — I 

For each MC generated configuration we first fix one of the 2 d flip sectors |l9|], where 
d are the spatial dimensions, and then fix Coulomb gauge maximizing separately for every 
time slice t the gauge functional: 

F 9 (*) = ^E Tr U i& ') ' U i& ') = 9& ') U i& £ ) 3\x + i, t) (3) 

with respect to local gauge transformations g(x,t) G SU(2), employing a simulated anneal- 
ing plus over- relaxation algorithm. The local maxima of F g (t) satisfy for each fixed t the 
differential lattice Coulomb gauge transversality condition for the gauge potentials: 

diA$(x,t) = Af(x + i/2,t) - Af(x -i/2,t) = 0. (4) 

We choose the best F g (t) out of n c random gauge copy for each time slice t separately, where 
usually n c = 5, and combine them into the best copy of the whole configuration for the 
flip sector chosen. We now go to the next sector and repeat the procedure. At the end we 
chose the sector with the highest global functional F g = ^2 t F g (t) as the best copy of the 
configuration. The number of total recursions is 2 d Ln c . Increasing n c leads of course to a 
better gauge fixing. However, going to higher volumes, one might consider to increase n c 
only for the sectors with the highest F g [19j. 

To extract the p dependence we need to fix the time gauge g(t). We choose here the 



integrated Polyakov gauge prescription proposed in Ref. 
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L 3 

corresponding in the continuum to 



u(t) = — U (x, t) — > const. (5) 



= y cPxdpAp&t) d J d 3 xA (x,t) = (6) 



To achieve this, we define 



u(t) 



u(t) 



e SU(2) 



(7) 



^Det [u(t)] 

Periodic boundary conditions in t make Tr]^[ 4 -u(t) = TrP invariant under g(t). Eq. (jSJ) can 
thus be fixed recursively through 



g(t)u(t)gi(t + l)=u = pV L 
up to a global gauge g(0) satisfying \q(0), P ] = 0. We can choose g(0) = t. 



(8) 



The main observation, made in Ref. 



21 



52[ for the 3+1 dimensional case and which we 



find here to hold in 2+1 dimensions as well, is that the bare gluon propagator factorizes as: 

M\P\) 9p{z) „ _ Po 



Dp(\p\,p ) 



(9) 



\p\ 2 1 + z 2 ' |p| 

The function gp(z) will in general depend on the temporal gauge g(t). However the static 
propagator, defined as 



urn 

ix\p\ 



dz 



l + z 



2 ; 



(10) 



will be independent of such choice. As shown in Ref. [21|, |52] lattice cutoffs to the integral 
in Eg . f[10l) cause scaling violations; as a viable 53j alternative to the lattice Hamiltonian 
limit |54j we define Dp(\p\) directly from fp(\p\), extracting it from the data after dividing 
out gp(z) and averaging over p : 



Dp(\p\) 



\p\ 



oc \p\ ^2Dp(\p\,po) 



l + z 2 



(ii) 



We obtain in this way a multiplicative renormalizable Dp(\p\), £) M (|p|) = Z(/3, p)Dp(\p\] 
We choose \i = oo, i.e. we fix the overall normalization such that limi^^oo |p|£) M (|p|) = 1. 
In 2+1 dimensions a Gribov like IR leading behavior: 



g(\p\: 



(12) 



plus corrections of the type: 



G(\p\) 



\P\ 



l+a 



(13) 



^/|^(2+a)/3 + yU (2+a)/3 

describes the data extremely well. To account for the intermediate momentum region differ- 
ent choices for the sub-leading terms other than Eq. (Tl3|) can be made, the only constraints 
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being asymptotic behaviors in the UV and o(\p\) in the IR. Lacking any theoreti- 

cal input further restricting the functional form is beyond our scope. In Fig. [I] we show 
Dc{\p\) = IpI^D^IpI), which nicely extrapolates to a constant in the IR, together with the 
fit to |p| _1 (cCr(|p|) + cG(\p\)). To improve readability versus the 3+1 dimensional propa- 
gator we have normalized the data to 3.5|p|~ 2 in the UV. We find M = 4.6(4) GeV and 
0-2 < x 2 /d.o.f. < 0.6, depending on the choice of the sub-leading terms, our main source of 
error. 

In 3+1 dimensions D^dpl) is well described by Gribov's formula alone 21] . In Fig. [TJwe 
show D c {\p\) = Ipl^D^lpl) and its fit to \p\- l G(\p\) as in Eq. flTJ); we obtain M = 0.856(8) 
GeV with % 2 /d.o.f. = 1.6. In Ref. {21] the deep IR data slightly deviated from a constant 
behavior, resulting in a xVd.o.f. = 3.3 for the same fit. These deviations are here, with 
increased volume and improved gauge fixing, considerably smaller, halving the x 2 /d-o.f. 
value. We therefore decided to ignore possible sub-leading corrections as in Eq. (fT3|) for 
the time being. Only a study on larger volumes and closer to the continuum limit, in the 



spirit of the Landau gauge analysis of Ref. [22, |23j, |55| , could clarify the situation. We also 



stress that no logarithmic corrections are necessary to describe the data in the UV. Although 
this is expected in 2+1 dimensions, it is a highly non trivial result in the 3+1 dimensional 
case, since one loop calculations for Dp(\p\,p ) give an explicit anomalous dimension [56 ]. 
However, contrary to covariant gauges, neither higher orders are under control nor a leading- 
log re-summation is viable. Cancellations when going from the one loop full propagator to 
the full static propagator are therefore in principle possible. 

As a last remark, both in 2+1 and 3+1 dimensions the asymptotic behaviors of D^dpl) 
are in agreement with the results of the IR and UV analysis of the DSE in the Hamiltonian 
approach of Ref. [15I \d7 \. 



B. Landau propagator 



The lattice gluon propagator in Landau gauge D L (p) has been the subject of extensive 
research; we refer the interested reader to the relevant literature regarding restoration of 
rotational invariance, gauge copy effects, continuum and thermodynamic limits etc. 16l-l20|. 



Here we simply apply the by now well established methods developed in the above references. 



In particular, we basically follow Ref. 
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IPI (GeV) 

FIG. 1. Dc(\p\) and fits in 2+1 and 3+1 dimensions. 

The results for 2+1 and 3+1 dimensions, normalized such that the peaks coincide with 
those of Fig. [TJ are shown in Fig. [2] and agree with all known previously published results. 

C. Landau vs. Coulomb 

On trivial dimensional grounds, based on their tree level UV behavior, it is obvious that 
D L {j>) should be compared to D c (\p\) = \~p\~ 1 D ^{\p\) . For both 2+1 and 3+1 dimensions 
a simple inspection of the Landau gauge transverse propagator plot against p = \/\p\ 2 + pi 
given in Fig. |2] shows striking similarities with Dc(\p\) plotted against \p\, as in Fig. CD 
The two curves almost coincide, up to a faster fall off of Di(p) in the UV. In fact, it 
simply looks as if the Landau 4-momentum scale p is stretched by some factor when com- 
pared to the Coulomb momentum scale \p\. For a quantitative comparison, motivated by 
this observation, we make the simple Ansatz D L {p) = D c (pl{p)), where Pl(p) = V pip/h) 
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2.5 




FIG. 2. Dl{p) in 2+1 and 3+1 dimensions. 



is strictly monotonous and p(p/A) is a dimensionless function interpolating between two 
constants pir and puv m 2+1 dimensions, where the anomalous dimension vanishes, and 
between a constant piR and puvlog 7//2 p 2 /AQCD in 3+1 dimensions. We use for p(p/A) a 
ratio of rational functions of the same order, p{p/A) = R n (p 2 / A 2 ) / Rd{p 2 / 'A 2 ) , corrected 
in 3+1 dimensions by log 1 ^ 2 Ri(p 2 / A q CD ), Ri being also a rational function. Polynomials 
of degree < 3 for R n ,d,i{x) are sufficient for our scope. In 3+1 dimensions shuffling the p 
dependence inside or outside the log is somewhat arbitrary. Indeed, a form of the type 
p(p/Aq CD ) = p IR log 7//2 Ri(p 2 /Aq CD ) is already sufficient, so that in both 2+1 and 3+1 di- 
mensions one single mass parameter, A or Aqcd is needed to describe p. A certain level 
of arbitrariness is also present in the normalization of the highest and lowest coefficient of 
Rn,d,i(%) as we ll as in their degree. The errors quoted below reflect all such, partly corre- 
lated, uncertainties, which will in particular show up in the value of A and Aqcd- However, 
lacking further theoretical input, pinning down p(p) to a precise functional form goes be- 
ll 



yond the scope of this work. We have fitted both |p|.Dc , (|p|) and Pl(p)Dl{pl{.p)) at the 
same time to Eq. ffT2|) (plus Eq. ffl3l) in 2+1 dimensions). p(p) as obtained from the fits is 
shown in Fig. [21 while the resulting comparison between D c (\p\) vs. D L {p) plotted against 
Pl{p) in 2+1 and 3+1 dimensions is given in Fig. H] and Fig. Numerically, we obtain 
p m = 1.44(14), p uv = 2.610(15), A = 3.0(3) GeV and xVd.o.f. = 1.2 in 2+1 dimensions, 
p m = 1.255(15), puv = 1-49(1), 7 = 0.52(3) and A QC d = 1.05(15) GeV with x 2 /d.o.f. = 1.6 
in 3+1 dimensions; the fitted values for all parameters in Eq. ffl2|) and Eq. ffl3|) . particularly 
for M, exactly coincide with those of Section III Al Moreover our 3+1 dimensional value for 
M = 0.856(8) GeV is in excellent agreement with recent SU(3) estimates of the IR mass in 



Landau gauge which explicitely include a dimensional two condensate in the analysis [58]. 

It is obvious that there is no way to relate Dc(\p\) to a scaling, IR vanishing Di(p) while 
also accounting for the UV region, since a monotonous mapping between the momenta 
could never do the job. One could of course object that the choice, in our opinion natural, 
to compare the Landau propagator with Dc{\p\) already shuts the doors to any scaling 
behavior. One could therefore attempt to compare Dl(p) directly to D^(\p\) again via a 
monotonous mapping of the momenta, e.g. -Dl(p) = D^p^p)), to check if an agreement to 
the same level of accuracy as in our analysis could be reached. In this case the function Pl{p) 
should behave up to dimensionful constants like oc p K in the IR and oc p 2 (plus logarithms 
in 3+1 dimensions) in the UV, i.e. Pl{p) would in principle need two new mass scales, 
besides the one for the intermediate momentum region, to account for the IR and UV 
behavior. As discussed above and as can be inferred from Fig. [21 this should be compared 
to the present analysis, where the only mass scale needed in p(p) basically coincides with 
the Gribov mass both in 2+1 and 3+1 dimensions, even though its exact value will depend 
on the, still arbitrary, parameterization of p(p). Moreover, that the Gribov parameter M 
extracted both from the lattice Coulomb and Landau gauge gluon transverse degrees of 
freedom coincide to such degree of precision constitutes in our opinion a highly non trivial 
result supporting a physical interpretation for the dynamically generated gluon mass. Of 
course the points raised above do not forbid in principle the existence of a mapping between 
the Gribov-like Coulomb propagator and a scaling like solution in Landau gauge, but lacking 
any theoretical arguments and as long as one can neither obtain a scaling solution on the 
lattice nor a solution obtained with continuum methods is available to sufficient degree of 
precision, any analogous quantitative check as the one performed in this paper to prove or 
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disprove such correspondence will remain impossible. 



T 




p (GeV) 

FIG. 3. p(p) in 2+1 and 3+1 dimensions, the latter compared to the perturbative asymptotic 
behavior. 

III. CONCLUSIONS 

We have shown that in 2+1 and 3+1 dimensional Yang-Mills theory, where physical gluons 
propagate, a natural and simple relation exists between the lattice results for the gluon-gluon 
Green's functions in two different gauges, i.e. the IR finite transverse Landau propagator 
Di(p) and the IR vanishing static transverse Coulomb propagator D^(\p\). Indeed, the two 
are equivalent if one simply rescales the momenta, i.e. Di(p) = Pl(p)^ 1 D i ^ 1 (pl(p)) where 
Pl{p) = P p{p) and p(p) is given in Fig.[3j D^dp]) agrees with the Gribov-Zwanziger scenario 
and has a natural interpretation within the dual-superconductor picture of confinement 
This corroborates a monopole condensation interpretation of the gluon massive solution 
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FIG. 4. Dc(\p\), Dl(j>l(p)) and fit in 2+1 dimensions. 



in Landau gauge as well, making, as long as no scaling like solution can be found on the 
lattice, the non perturbative loss of BRST symmetry more bearable. Whether the static 
Coulomb gluon itself coincides with the Yang-Mills field's physical degrees of freedom is 
still an open issue 



59 



601 ] . Should the uniqueness of the IR massive gluon behavior on the 



lattice be confirmed, it would suggests that in covariant gauges in more than one spatial 



dimension the Gribov-Zwanziger 



3| term should be explicitly included in the action, 



BRST symmetry would be non-perturbatively broken and the Kugo-Ojima criterion 



could not be fulfilled. Indeed an effectively massive gluon propagator is natural in a theory 
where monopoles condense, as can be explicitly seen in the Georgi-Glashow model |6ll 
621 ] : dual-superconductivity is expected to work in a similar way to confine gluons, albeit 
with a dynamically generated scalar boson |g3 | . Indeed, a topological origin for dynamic 
mass generation in QCD has been argued long time ago in Ref. 64J and offers a natural 



interpretation of the pinch technique result (9 



101 ] . In topological theories on the other 
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2 4 6 8 10 12 

Ipl. P L (GeV) 

FIG. 5. Dc(\p\), Dl{pl(p)) and fit in 3+1 dimensions. 

hand, where the physical Coulomb propagator should identically vanishes, a IR conformal 
confinement with a non perturbatively conserved BRST charge can be realized. The Chern- 
Simons and BF theories will be interesting cross-checks to confirm the different mechanism 
between topological and QCD like theories. 
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